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Stochastic  Models  for  the  Random 
Location  of  Individuals  in  a Habitat 


by 

Naftali  A.  Langberg,  Roger  Johnson,  and  Ralph  A.  Bradley 

ABSTRACT 

Very  little  had  been  done  in  the  past  to  develop  coherent  stochastic 
structures  for  the  location  of  individuals  in  a region.  Such  a structure 
can  provide  sound  probabilistic  bases  for  a variety  of  inference  procedures 
presented  in  the  literature  and  clear  the  way  to  the  development  of  consis- 
tent statistical  methodology.  Various  probabilistic  models  for  the  random 
locations  of  individuals  in  a region  are  considered. 

Models  for  both  infinite  and  finite  regions  are  developed.  In  the  liter- 
ature, certain  marginal  distributions  are  used  to  characterize  the  numbers 
of  individuals  in  subregions.  Relationships  between  these  marginal  distribu- 
tions and  our  probabilistic  models  are  discussed. 


Key  words:  Stochastic  processes,  stationarity,  independent  increments,  Poisson 
Process,  Poisson,  Binomial  and  Negative  Binomial  distributions. 
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1.  Introduction  and  Summary. 

Trees  in  a forest,  stars  in  a galaxy  and  bacteria  on  a Petridish  are  a 
few  examples  wherein  there  is  interest  in  the  spatial  distributions  of  indi- 
viduals over  regions.  In  pursuit  of  better  understanding  of  these  distribu- 
tions, attempts  have  been  made  to  use  probabilistic  and  statistical  tools 
[Eberhardt  (1967),  Pilou  (1969),  Pollard  (1971),  Patil  and  Stiteler  (1974), 

Cox  (1976),  Cox  and  Lewis  (1976)].  Very  little  effort  has  been  made  to  develop 
coherent  statistical  structures  for  the  locations  of  individuals  in  a region. 
The  main  purpose  of  this  paper  is  to  propose  such  structures  in  order  to 
avoid  inconsistencies  that  can  be  found  in  the  literature,  to  provide  a sound 
probabilistic  basis  for  a variety  of  inference  procedures  that  are  used,  and 
to  permit  the  development  of  better  and  consistent  statistical  methodology  in 
the  area. 

In  Sections  2,  3 and  4,  we  show  that,  under  certain  conditions,  the 
collection  of  random  variables  representing  the  numbers  of  individuals  in 
subsets  of  a region,  can  be  regarded  as  a stochastic  process.  All  the  models 
for  the  stochastic  process  presented  here  have  a stationarity  property  defined 
in  (2.1)  below.  In  Section  2 we  present  Simple  Random  Models.  The  distribu- 
tion of  the  number  of  individuals  in  a region  is  considered  to  be  simple  random 
if,  in  addition  to  stationarity,  it  has  a property  of  independent  increments 
defined  in  (2.2)  below.  We  prove  that  simple  random  models  are  uniquely  deter- 
mined by  a positive  number  and  a sequence  of  real  numbers  in  [0,  1]  summing 
to  unity.  In  a number  of  the  references  cited  above,  Poisson  and  Negative 
Binomial  distributions  have  been  used  to  describe  the  number  of  individuals  in 
subsets  of  a region.  We  conclude  Section  2 by  showing  that  these  two  distribu- 
tions are  marginal  distributions  resulting  from  particular  siiiq>le  random  models. 
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In  Section  3 we  introduce  models  that  are  stationary,  but  do  not  necessarily 
have  independent  increments.  We  show  that  a Negative  Binomial  distribution 
proposed  by  Patil  and  Stiteler  (1974)  to  characterize  the  random  behavior  of 
the  numbers  of  individuals  in  subsets  of  a region  is  a marginal  distribution 
resulting  from  two  different  stochastic  models  that  do  not  have  independent 
increments. 

Section  4 is  devoted  to  the  development  of  several  different  stochastic 
models  when  the  region  under  consideration  has  finite  Lebesgue  measure.  It 
is  shown  that  the  Binomial  distribution  suggested  by  Eberhardt  (1967)  can  be 
derived  from  of  these  models. 
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2.  Simple  Random  Models. 

Some  basic  notation  is  needed  to  define  spatial  stochastic  models.  Through- 
out the  paper,  R will  denote  a Euclidian  space  or  a Borel  subset  of  a Euclidian 
space  with  finite  Lebesgue  measure,  A will  be  the  Lebesgue  measure  on  R,  X(A) 
will  be  the  number  of  individuals  in  a subset  A e R,  and  ^ will  represent  the 
set  of  all  Borel  subsets  of  R with  finite  Lebesgue  measures. 

Individuals  are  said  to  be  located  in  R in  a simple  random  way  if  defined 
propositions  (2.1)  and  (2.2)  hold. 

StationArity;  For  every  set  A e the  distribution  of  the  random  variable  X(A) 

(2.1) 

depends  only  on  A (A) . 

Independent  increments;  If  A. , ...,  A represent  any  m disjoint  sets  in  J, 

(2.2) 

m > 1,  the  random  variables,  X(A,),  ....  X(A_)  are  Independent. 

The  objective  of  this  section  is  to  show  that  under  conditions  (2.1)  and  (2.2), 
the  collection,  ^ stochastic  process  determined  by  a positive  num- 

ber and  a sequence  of  nissbers  in  [0,  1]  summing  to  unity. 

Let  q be  an  arbitrary  point  in  R,  Y^(t),  the  number  of  individuals  in  a 
sphere  centered  at  q having  Lebesgue  measure  t,  t i 0,  Y^(0)  = 0,  A^,  A2,  sets 
in  and  k^,  k2  two  nonnegative  integers.  Then 


P(X(Aj)  - kj,  X(/.2)  * 1-2: 


I P[X(A  n A2)  - j.  X(A  n = r,  X(A  n 
j*r-k 
j+t-kj 


1].(2.3) 


From  (2.1),  (2.2)  and  (2.3)  it  follows  that 


P[X(Aj)  - kj,  X(A2) 


k2l 


jT-k  **  ‘ 


j]  P[Y  (t2)  - r]  P[Y  (tj)  - t]. 


j*i«k: 


(2.4) 
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where  t2  and  are  X(Aj  n A2),  X(Aj  n A2)  and  X(Aj  n A2)  respectively. 
Without  difficulty,  equations  (2.3)  and  (2.4)  can  be  extended  to  any  finite 
number  of  sets  in  The  formal  proof  is  omitted  but  the  general  result  is 
stated  in  the  following  theorem. 


Theorem  2.1.  Under  conditions  (2.1)  and  (2.2),  the  joint  distribution  of  any 
finite  collection  of  random  variables  from  tX(A)}^^^  is  determined  uniquely  by 
a joint  distribution  of  a finite  collection  of  random  variables  from 


To  achieve  our  objective,  it  suffices  to  show  that,  under  (2.1)  and  (2.2),  the 
collection  is  a stochastic  process  determined  by  a positive  number 

and  a probabilistic  sequence. 

Conditions  (2.1)  and  (2.2)  respectively  insure  that  ^YqCt))^^Q  is 
stationary  and  has  independent  increments.  We  state  formally: 

For  every  two  nonnegative  numbers  t and  s,  Y^(f)  and  Yq(t  s)  - Yq(s)  are 
identically  distributed.  For  every  m real  numbers  0 5 tj  s t2  s ...  s t^^^, 
the  random  variables  Y^(tj),  Y^(t2)  “ inde- 

pendent. Applying  a well  known  characterization  of  stationary  processes  with 
independent  increments,  [Khintchine  (1960),  36.],  we  obtain  the  following: 

Theorem  2.2.  Under  conditions  (2.1)  and  (2.2),  there  exist  a Poisson  process 

{N(t))^^Q  and  a sequence  of  i.i.d.  random  variables  2 

positive  integers,  independent  of  {N(t)}  such  that  the  processes  {Y_(t)} 

t2o  q tao 

and  { > Z, )..-  are  identical. 
l‘l  ItiO 
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It  is  a consequence  of  Theorem  2.2  that,  the  collection  * 

stochastic  process,  determined  by  the  constant  E[N(1)]  and  the  sequence 


(P(Z,  - i)> 


i-1,2,...* 


Two  marginal  distributions  for  X(A),  used  frequently  in  the  literature, 
arise  from  special  cases  of  simple  random  models;  they  are  the  Poisson  and 
the  Negative  Binomial.  To  show  it,  we  give  the  two  definitions: 


X has  a Negative  Binomial  distribution  with  parameters  a and  b,  (a,  b > 0)  if 

la 


(r^)  (rrr)’-  ^ ^ “»•«  >•(«) 

m 

/e”*x°"^dx,  o > 0. 

0 


X has  a logarithmic  distribution  with  parameter  0 (0  < 9 < 1) , if  P(X  ■ i)  ■ 
-eViinCl  - 0),  i « 1,  2,  ...  . 

If  Zj  of  Theorem  2.2  is  a degenerate  random  variable  at  1,  or  has  the  logarithmic 
distribution  defined  above,  it  follows  from  Theoran  2.2  that,  for  every  set  A e 
the  distribution  of  X(A)  is  respectively  Poisson,  with  mean  X(A)  E[N(1)],  or 
Negative  Binomial  with  parameters  a ■ ■ tn(P-'^9)  ' ^ V e 
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3.  A Generalization  of  the  Simple  Random  Models. 

In  this  section  we  present  two  methods  of  developing  a general  stationary 

t ^ ■ ■ 

stochastic  structure  for  the  random  quantities  This  is  dppe  in  twi 

steps.  Firstly,  we  present  joint  distributions  of  X(Aj),  X(A^),  where 
Aj,  A^  are  any  m disjoint  sets  of  I and  m = 1,  2,  ...  . We  then  show 
that  this  set  of  finite  dimensional  distributions  determine  uniquely  the 
stochastic  structure  of 

For  construction  purposes  only  we  introduce  some  notation. 
be  a stochastic  process  on  the  nonnegative  integers,  stochastically  increasing 
in  t,  that  is,  M(t)  s M(s)  for  0 s t s s and  {W. _ a sequence  of  random 

X X ^ ^ • • • 

variables  with  positive  integer  values,  independent  of  We  assume 

that  and  2 depend  on  parameters  0j^  and  02  respectively, 

ranging  in  parameter  spaces  and  q^,  where  ©2  is  ^ collection  of  infinite 
se",uences,  0j  and  ©2  being  probability  spaces  with  F and  G the  respective 
probability  measures. 

Let  Aj,  ...,  A^  be  disjoint  sets  in  ^ with  respective  Lebesgue  measures 
^1'  *”*  ^m  *'l’  **'’  ^m  nonnegative  integers.  We  define  the  joint 

probability  in  two  ways. 

m Mftj) 

P[X(A  ) - k , i - 1,  ...,  m]  = / / n P[  J W = kJdF(0  )dG(02)  (3.1) 

®l’‘®2^“^  j=l  ^ 

or 

Mfti) 

P[X(A  ) - k.,  i - 1,  ...,  m]  - / / P[  ^ W = k.,  i = 1,  ...,  m]dF(0,)dG(0,).  (3.2) 

' ^ 0lX02  j=l  ^ ^ 1 Z 
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The  consistency  of  the  set  of  joint  distributions  of  finite  collections  of 
random  variables  from  generated  by  disjoint  sets  and  defined  by 

(3.1)  or  (3.2),  follows  from  the  stochastic  structure  imposed  on 
< {tl(t)}.^n  {W.}.  , „ The  stationarity  condition,  stated  in  (2.1),  is 

clearly  satisfied,  since  definitions  (3.1)  and  (3.2)  vary  only  with  the 
Lebesgue  measures  of  the  respective  sets.  The  generalization  of  equation  (2.3) 
to  any  finite  number  of  sets  in  ^ provides  a way  to  extend  definitions  (3.1) 
and  (3.2)  to  any  finite  number  of  sets  in  J.  These  extensions  preserve  the 
required  consistency  and  stationarity  of  ^X(A)}^^p  We  have  proved  the 
following: 

Theorem  3.1.  If  the  joint  distributions  of  every  finite  number  of  random 
variables  in  tX(A)}^^^  are  given  by  the  extensions  of  (3.1)  or  (3.2),  then 
{X(A)}^^j  is  a stationary  stochastic  process. 

It  is  important  to  note  that  (3.1)  and  (3.2)  yield  the  same  marginal  distribu- 
tions for  X(A),  A e but  they  define  different  processes.  Statistical 
methodologies  based  on  the  two  models  may  be  quite  different  and  this  has  been 
ignored  in  the  literature. 

In  Section  2,  we  proved  that  simple  random  models  depend  on  a positive 
number  say  w,  and  a probabilistic  sequence,  say  {P,}.  , _ . Now  we  show  that 

those  models  are  particular  cases  of  the  stochastic  processes  presented  by  (3.2). 

To  do  so,  let  {M(t)}  « be  a Poisson  process  with  parameter  y,  and  (W. }.  , - 
be  an  i.i.d.  sequence  of  random  variables  given  by  P(Wj  * i)  = Pj^,  i > 1,  2,  ...  . 

In  addition  we  take  3^  and  02  to  be  sets  containing  only  y and  the  sequence  2 

respectively.  Now  the  extension  of  (3.2)  reduces  directly  to  a single  random 


model. 
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Let's  assume  that  the  sequence  of  discrete  random  variables  {W.}.  , - 
degenerate  at  1,  0^  is  [0,  “) , ©2  is  only  the  sequence  (1,  0,  0,  ...)» 
{M(t)}^^Q  is  Poisson,  and  F is  given  by 


IS 


F(x)  = 


_1_  f v“"^dv 

r(o)  J ® y 


X > 0 


X < 0 


a,  8 > 0. 


(3.3) 


Then  the  marginal  distribution  of  X(A),  A e computed  by  (3.1)  or  (3.2),  is 
negative  Binomial  with  a = a and  b --  A (A).  This  marginal  distribution  for 

P 

X(A)  was  assumed  by  Eberhardt  (1967)  and  by  Patil  and  Stiteler  (1974). 
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4.  Models  for  finite  regions. 

In  the  previous  two  sections,  we  presented  stochastic  models  when  indi- 
viduals were  randomly  located  in  an  infinite  region.  In  this  section,  five 
different  models  are  developed  for  situations  wherein  individuals  are  located 
in  a Borel  subset  of  a finite  Lebesgue  measure  in  R. 

Let  be  the  set  of  all  Borel  sets  contained  in  R^  and  X^(A),  the 
number  of  individuals  in  set  A,  A e Our  objective  is  to  develop  a 

stochastic  structure  for  the  collection  • One  way  to  acheive  this 

objective  is  to  select  a stochastic  process  developed  in  Section  3 for  the 
infinite  region  R,  say  condition  it  by  an  event  related  to  the 

random  variable  X(R^).  To  be  move  formal,  let  A^,  . . . , A^  be  sets  in  and 
let  kj,  ...,  be  nonnegative  integers.  We  define  the  desired  probability 

P[Xf(Ap  = k^,  i = 1,  ...,  m]  = P[X(Ap  = k.,  i = 1,  ...,  m|X(R^)].  (4.1) 

The  consistency  and  stationarity  of  the  process  ^X^(A)}^^^  defined  in  (4.1) 
are  self-evident. 

For  the  four  remaining  alternative  ways  of  developing  a stochastic 
structure  for  ^X£.(A)}^^^  , we  refer  to  the  process  {^(t)}^^^,  the  sequence 
{W.}.  , , , the  sets  ©,  and  and  the  probability  measures  F and  G 

which  were  introduced  in  Section  3.  We  define  first  the  joint  distributions 
of  X£(Aj),  ...,  X^(A^),  where  A^^,  ...,  A^  are  disjoint  sets  in  and 
m » 1,  2,  ...,  then  extend  these  definitions  to  any  finite  collection  of 
sets  in  Since  the  extension  technique  has  been  used  twice,  it  is  omitted 

from  our  present  discussion.  Now  let  A^,  . . . , A^  be  disjoint  subsets  of  R^, 
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= tp,  A(A^)  = t^,  i = 1,  m,  and  let  k^,  ...,  k^  be  nonnegative 
integers.  We  define  the  desired  probability  in  the  following  four  ways: 

m mCt-) 

P[X-(A  ) = k . i = 1,  m]  = / / n P(  2 W = k I ^ W.)dF(0,)dG(0,). 

ejxe2i=i  j=l  J j=l  J ^ ^ 

P[X£(A.)  = k , i = 1,  ....  m]  = / f P(  I W = k. , i = 1,  m|  2 W.)dF(e, 

01X02  j=l  ^ ^ j=l  J 1 


(4.2) 


)dG(02). 

(4.3) 


m Mftj) 

m]  = / / n P(  [ W = k |M(t_))dF(0  )dGC0-). 

0^x0  i=i  j = i J 1 A ^ 


P[Xf(Ap  = k^,  i = 1,  ....  mj  = j j H PC  I W.  = kJM(tJ)dF(0,)dGC0,).  (4.4) 

M(ti) 

P[Xf(A.)  = k , i = 1,  ....  m]  = / / P(  [ W = k. , i = 1,  ....  m|M(t  'dF(0.)dG(0  ). 

9,x0  j = i J A 0 1^ 

(4.5) 


The  stochastic  processes  determined  by  (4. 2) -(4. 5)  have  been  constructed  to 
be  both  consistent  and  stationary.  Definitions  (4.2)  and  (4.3)  or  (4.4)  and 
(4.5),  yeild  identical  marginal  distributions  for  X^(A),  A e but  define 
different  processes. 

If  we  let  (M(t)}^>Q  be  a Poisson  process,  (W.}.  _ , , , a sequence 

of  random  variables  degenerate  at  1,  0^  and  0^  singleton  sets,  then  the 
marginal  distribution  of  X^(A),  A e according  to  each  of  the  last  four 
definitions  are  Binomial  with  X(R£)  corresponding  to  the  number  of  Bemulli 
trials  and  corresponding  to  the  probability  of  succes  in  a single  trial. 
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5.  Concluding  Remarks. 

The  selection  of  models  for  stochastic  processes  for  particular  appli- 
cations and  the  development  of  pertinent  statistical  methodologies  have  not 
been  addressed  in  this  paper.  V/e  have  demonstrated  that  models  proposed  may 
be  used  to  yield  marginal  distributions  assumed  in  the  literature;  they  may 
also  be  used  to  avoid  unwarranted  assumptions  and  inconsistencies  that  arise. 
We  propose  in  subsequent  work  to  use  the  general  models  of  this  paper  to 
devise  improved  statistical  methodologies  for  problems  involving  the  location 
of  individuals  in  a habitat. 
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development  of  consistent  statistical  methodology.  Various  probabilistic  models  for  the 
random  locations  of  individuals  in  a region  are  considered. 

Models  for  both  infinite  and  finite  regions  are  developed.  In  the  literature,  certain 
marginal  distributions  are  used  to  characterize  the  numbers  of  individuals  in  subregions. 
Relationships  between  these  marginal  distributions  and  our  probabilistic  models  are 


discussed. 


